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Abstract 

We show that smoothing of multiaffine surfaces that are generated by 
simulating a crosslinked polymer gel by a frustrated, triangular network of 
springs of random equilibrium lengths [G.M. Buendia, S.J. Mitchell, P.A. 
Rikvold, Phys. Rev. E 66 (2002) 046119] changes the scaling behavior of 
the surfaces such that they become self-affine. The self-affine behavior is 
consistent with recent atomic force microscopy (AFM) studies of the surface 
structure of crosslinked polymer gels into which voids are introduced through 
templating by surfactant micelles [M. Chakrapani, S.J. Mitchell, D.H. Van 
Winkle, P.A. Rikvold, J. Colloid Interface Sci., in press]. The smoothing 
process mimics the effect of the AFM tip that tends to flatten the soft gel 
surfaces. Both the experimental and the simulated surfaces have a non-trivial 
scaling behavior on small length scales, with a crossover to scale-independent 
behavior on large scales. 

PACS: 61.43.Hv; 89.75.Da; 82.70. Gg; 68.37.Ps 
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1 Introduction 

In a recent study jl] we proposed a frustrated spring-network model to simulate 
the surface structures of crosslinked polymer gels into which voids are introduced 
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through templating by surfactant micelles. Experimentally, templated polyacry- 
lamide gels are created by mixing acrylamide with a crosslinker in the presence of a 
surfactant p]. The resultant gels have a wide range of pore sizes, making them ideal 
materials for macromolecular separation [Hj. The crosslinking process introduces 
inhomogeneities in the spatial density of the gel that affect its surface configuration. 
Our results show that the surfaces generated by the spring-network model have a 
nontrivial, multiaffine scaling behavior at small length scales, with a crossover to 
scale-independent behavior on large scales. 

At small length scales we found that the gth order increment correlation func- 
tion has a power-law behavior, with g-dependent exponents that change with the 
vacancy concentration. However, experimental surfaces of surfactant-templated hy- 
drogels studied by atomic force microscopy (AFM) 4J show g-independent self-affine 
scaling at length scales below the crossover length. In Ref. we discussed the pos- 
sibility that the multiaffine behavior was due to a finite density of vertical disconti- 
nuities due to overhangs in the simulated surfaces. Later work showed analytically 
and numerically that when vertical discontinuities are introduced into a self-affine 
surface, it becomes multiaffine Here we extend this study to hydrogel surfaces, 
and we show that smoothing the simulated surfaces by convoluting them with a 
Gaussian changes the scaling behavior from multiaffine to self-affine, in agreement 
with the experimental AFM results. The smoothing process mimics the effect of the 
measurement process of an AFM tip interacting with the soft substrate. 

2 The frustrated spring-network model 

The model, described in detail in Ref. P, consists of a two-dimensional triangular 
network of nodes interconnected by massless harmonic springs. The nodes represent 
the crosslinkers, and the springs correspond to the polymer chains. We impose 
periodic boundary conditions in the horizontal direction, while the top nodes are 
free to move in two dimensions. The bottom nodes are fixed, corresponding to 
bonding to a rigid substrate. 

We assume that each spring corresponds to a collapsed polymer chain with equi- 
librium end-to-end length /oi, which is randomly chosen with probability density 



Thus, loi is proportional to the square root of an exponentially distributed number of 
monomers between crosslinkers jSj. The constant 7 is proportional to the inverse of 
the average number of monomers between crosslinkers. Consistent with the picture 
of a random-coil collapsed polymer chain, we require that the elastic constant of 
each spring, fcj, be proportional to its equilibrium length, fcj = l/loi in dimensionless 



For each statistical realization of the system, the equilibrium spring lengths are 
chosen independently according to Eq. (1). In the initial network configuration, all 
springs are placed on the bonds of a regular triangular lattice of unit lattice constant, 
such that all springs are stretched or compressed to unit length, h = 1. This initial 





units |7j. 



2 



configuration is highly stressed. However, by choosing the value of 7 such that the 
average force exerted by each spring is zero, we ensure that the initial configuration 
is not globally stressed. From the initial configuration, the system is relaxed to one 
of its many local energy minima by a limited-memory Broyden-Fletcher-Goldfarb- 
Shanno (L-BFGS) quasi- Newton minimization algorithm |H|, and by repeating the 
calculations for several different realizations of the equilibrium spring lengths we 
verified that the statistical properties of the different relaxed configurations are 
essentially the same. 

Spring networks of x Ly = 1024 x 768 nodes were generated, and voids 
were introduced by randomly removing nodes and the springs connected to them. 
The Hoshen-Kopelman algorithm jH] was applied to eliminate clusters that were 
not connected to the fixed substrate, thus mimicking the removal of non-bonded 
material which is washed away during the experimental preparation. 

3 Calculations 

The one-dimensional surfaces are defined as the set of surface heights, h{xi), at 
equally spaced discrete horizontal points, Xi. The surface height at Xi, h{xi), is taken 
as the largest vertical distance, measured from the bottom, of the intersections of 
the springs with a vertical line at Xi, where each spring is represented by a straight 
line segment between the connected nodes. This definition of the surface is described 
in more detail in Ref. [Tj. 

The gth order generalized increment correlation function is defined as 

C,{r) = {\hixo + r)-h{xo)n (2) 

For many surfaces, Cg displays a power-law behavior, such that 

Cg{r) oc r'^"' (3) 

where Hg is the generalized Hurst exponent fOJ ■ For self-affine surfaces, the 
scaling behavior is independent of q, Hg = H, but for multiaffine surfaces, there is 
an infinite hierarchy of scaling exponents, such that Hg depends continuously on q, 
at least for some range of q values and length scales jT2I- 

Very often surfaces present a mixed behavior. For length scales smaller than some 
crossover, r < l^, Cq displays a non-trivial scaling behavior, while at length scales 
larger than this crossover, r > l^, Cq is constant, Cg{x) ~ Cg'^^. In our previous 
work [1^ it was shown that the simulated gel surfaces are multiaffine, and the scaling 
exponents, crossover length scales, and saturation values of the increment correlation 
function are strongly dependent on the void concentration. Despite the fact that the 
experimental surfaces seem to have a similar dependence on the void concentration, 
one major difference between the simulated and experimental surfaces remains: The 
experimentally measured surfaces are clearly self-affine on short length scales [21 E] , 
while the simulated surfaces are clearly multiaffine P?^ 

Recent work |5j indicates that vertical discontinuities are one important source of 
multiaffinity, and the surfaces generated by our model certainly have many vertical 
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Figure 1: A typical portion of the spring-network surface (20% vacancies) before 
and after smoothing by convolution with a Gaussian. The smoothing removes the 
vertical discontinuities without significantly altering the overall surface features. 
Such smoothing should mimic the measurement process of using a weakly interacting 
AFM tip to probe the structure of a soft surface. 



discontinuities caused by overhangs. Even with a weakly interacting AFM tip, 
experimental measurement of a soft surface, such as the polyacrylamide gels studied 
in should have the effect of smoothing out any vertical discontinuities and 
overhangs, thus showing a self-affine surface, even if the real surface is multiaffine. 
To mimic the effect of passing a weakly interacting AFM tip over the spring-network 
surfaces, we have smoothed the simulated surfaces by convolution with a Gaussian 
of cr = 0.1 in the dimensionless spatial units of the problem. This process eliminates 
the overhangs without changing the overall structural properties of the surfaces. 

Figure 1 shows a detail of the generated surface and the smoothed surface. Fig- 
ure 2 shows the gth root of the generalized increment correlation function Cg^'^{r) 
for systems with 0%,40%, and 49.5% vacancies, for both the original and smoothed 
surfaces. It is evident that when the discontinuities are eliminated by smoothing, 
the Hq values become independent of q, thus indicating that the surfaces are no 
longer multiaffine but have become self-affine after smoothing. The Hg exponents 
are plotted in Fig. 3 versus the void concentration for the original spring-network 
surfaces. For the smoothed surfaces, Hg = 1.0 for all void concentrations, as ex- 
pected for a smooth surface. The multiaffine scaling observed is consistent with that 
caused by vertical discontinuities [5j. 

4 Conclusions 

In this paper we have shown that surfaces generated by a frustrated spring-network 
model can reproduce several aspects of the scaling behavior of surfaces of surfactant- 
templated polyacrylamide gels, as observed in recent AFM experiments 0. Both 
real and simulated surfaces show a nontrivial scaling behavior characterized by a 
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Figure 2: The gth root of Cq(r) from the spring-network surfaces for q = 0.5 to 
q = 4.0 in steps of 0.5. For graphical simphcity, the q labels have been omitted, but 
Cj(r) > Cj{r), when i > j. The surfaces are averaged over 8, 8, and 10 independent 
reahzations of the equihbrium spring lengths for the networks with 0%, 40%, and 
49% vacancies, respectively. The linear regions of the plot indicate power-law scaling, 
and for the original simulated surfaces, the different power-law exponents (slopes in 
the log-log plots) indicate g-dependent multiaffine scaling. After smoothing, the 
power-law scaling at small length scales is self-affine (parallel lines in the log-log 
plots) with Hq — 1 for all q. 
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Figure 3: Hg vs. 1/q for the spring-network surfaces. Open symbols indicate 
the original, multiaffine surfaces, while corresponding filled symbols indicate the 
smoothed, self-affine surfaces. As expected, Hg = 1.0 for all of the smoothed sur- 
faces. The observed multiaffine scaling of the original surfaces is consistent with the 
behavior caused by vertical discontinuities as expected from Ref. 0. The data are 
averaged over 8 independent realizations of the equilibrium spring lengths for 0% 
to 47% vacancies, 10 realizations for 49% vacancies, and 14 realizations for 49.5% 
vacancies. 

power-law form of the increment correlation function at small and intermediate 
length scales. At larger length scales the increment correlation function reaches a 
saturation value. This large-scale behavior reflects the average network structure, 
while the small-scale behavior is due to microscopic density fluctuations. Below 
some crossover length, the experimental surfaces are self-affine while the original 
simulated surfaces are multiaffine. We show that after smoothing by convolution 
with a narrow Gaussian, which eliminates vertical discontinuities, the simulated 
surfaces become self-affine. Such smoothing mimics the effect of passing a weakly 
interacting AFM tip over the soft hydrogel surfaces. The large-scale behavior is not 
affected by the smoothing process. 
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